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I. INTRODUCTION 





A number of new tasks undertaken during the semi- 
annual reporting period are discussed (Sections 2 and 3) 
A final report on the H 2 -kinetics modeling work is 
contained herein as Section 4. Our continuing work 
dealing with turbulent reaction rate modeling is in 

Section 5. An addition to these items, two publications 
are in preparation: 

1. Chinitz, W. and Evans, j. s., "a Model for 

Reaction Rates in Turbulent Reacting Flows”, NASA Tech. 
Memorandum. 

2. Chinitz, W. , "Models for Calculating Reaction 
Rates in Turbulent Reacting Flows", paper to be presented 
at the 1984 JANNAF Propulsion Meeting, Feb. 7-9, 1984, 

New Orleans, LA. 




2. ALTERNATIVE SILANE COMBUSTION MECHANISM 



The use of silane CSiH^) as an effective ignitor and 
flame stabilizing pilot fuel has been well documented (e.g. ref. 1) . 
In order to better understand its behavior at the conditions en- 
countered in the combustor of a SCRAMJET engine, and to be able 
to predict its behavior at conditions which cannot be readily 
achieved in the laboratory, a reliable chemical kinetic mechanism 
must be developed. The first attempt at the formulation of such 
a mechanism for SiH^ is in ref. 1. It was shown subsequently 
(ref. 2), however, that the mechanism in ref. 1 yields predicted 
ignition delay times which are about five times longer than those 
experimentally obtained in a shock tube. 

As a consequence, an alternative mechanism was developed 
in ref. 3 so as to obtain agreement with the experimental shock 
tube data. The mechanism in ref. 3 postulates that silane pyrol- 
ysis 

SiH. -*• SiH- + H- ... 

422 (l) 

along with its direct oxidation 

SiH 4 + 0 2 -* SiH 3 + H0 2 (ii) 

are the initiating reactions. An examination of the magnitudes 
of the rate constants recommended in ref. 3 indicates that 
reaction (i) , along with the extremely rapid subsequent oxidation 
of Si» 2 

SiH 2 + 0 2 - HSiO + OH, (iii) 

are principally responsible for the disappearance of silane at 
SCRAMJET engine operation conditions. • . 



An alternative mechanism has been developed wherein the 
reaction 

SiH 4 + 0 2 SiH 2 0 + H 2 0 (iv) 

along with the hydrogen atom abstraction reaction 

SiH 4 + M SiH 3 + H (v) 

are the principal initiating reactions. Reaction (iv) is 
likely to be semi-global in nature, similar to the serai-global 
reactions in the methane mechanism in ref. 4. The complete 
silane oxidation mechanism postulated is in Table I. Complete 
details of its formulation will be presented in a forthcoming 
NASA Contractors Report. In this status report, we focus upon 
the results obtained using this mechanism. 

A comparison between the shock tube data and the theoreti- 
cal ignition delay times calculated using the proposed mechanism 
is in figs. 1 and 2. Generally, excellent agreement is obtained 
over the range of experimental conditions. A comparison between 
this mechanism and that in ref. 3 obtained for a number of 
SiH 4 /H 2 mixtures at one atmosphere, equivalence ratios of 1.0 
and 0.5, in the temperature range from 800K to 1250K, is in figs. 
3 and 4. For most of the conditions examined, the proposed 
mechanism predicts substantially lower ignition delays. Only for 
the 2% SiH 4 mixture at initial temperatures exceeding about 930K 
does the mechanism in ref. 3 predict shorter ignition delays. 

The predicted effect of equivalence ratio on ignition delay 

< 6 < 

time is in fig. 5. As can be seen, in the range 0.5 — — 10.0, 

increasing results in reduced ignition delays. This may parti- 
ally explain the empirically-observed pyroforicity of silane at 
one atmosphere, room temperature conditions since, under the 
circumstances in which this phenomenon has been observed, highly 



fuel-rich regions are certain to exist. The effect of decreas- 
ing the pressure to 0.5 atm is shown in fig. 6. As would be 
anticipated from prior chemical kinetics studies, decreasing the 
pressure results in increased delay times for SiH^/Hj mixtures 
down to 2% silane. 

Perfectly stirred reactor calculations were carried out to 
assess the flame-stabilizing behavior of silane. In fig. 7, 
blowout limit correlations are shown for pure silane at 1.0 atm, 
600K initial temperature using the proposed mechanism and that 
in ref. 3. As can be seen, the mechanism in Table I predicts a 
reduced stable flame region when compared with that in ref. 3. 
Both, however, exhibit the same general shape, with (m/V)max 
occurring at about <J> =2.5, rather than at the anticipated value 

of <f> = 1.0. The reasons for this behavior, and whether or not 
it is in fact real, remain a matter of speculation. Experimental 
PSR studies would be required to help settle these matters. 

Computed blowout limit correlations using the mechanism in 
Table I are shown in figs. 8-10 for 1.0 atm, T q = 300K, 600K and 
1200K, respectively, for 100% SiH^, 20% SiH^/80% Hj and 100% • 

The ifigures make clear that the principal benefits to be derived 
from using silane for purposes of flame stabilization are derived 
using larger silane concentrations in regions of high equivalence 
ratio. In low $ regions, pure hydrogen produces a larger stable 
flame region. Also of interest is the unusual "reversal" behavior 
at T q = 1200K in which, for the 20% SiH^/80% mixture, for 
example, above 4> ** 1.0, increasing equivalence ratio permits 
higher values of m/V. Similar behavior is observed for 100% 




SiH 4 down to $ « 0.5. The reality of this behavior must also 
await experimental PSR results. 

The results in figs 8-10 are shown once again in figs. 

11-13 as $ versus residence time. The effect of pressure on 
blowout limit behavior is shown in figs. 14 and 15. As expected, 
reducing the pressure tends to reduce the stable flame region. 

Consideration of computer run times and storage limitations 
require the development of a global silane oxidation model. 

Such a model has been formulated for the mechanism in Table I. 
The model consists of the three reactions: 

2 SiHj + 0 2 = 2 HSiO + 3H 2 

2 HSiO + 0 2 = 2 Si0 2 + H 2 

2 H 2 + 0 2 = 2 H 2 0 

with the rate coefficients given by 


(vi) 

(vii) 
(viii) 


(0.95234 <f> + - - 4 ] 62 - - 4.4286) 


‘f ,vi 


5.644 


1015 T 2 - 30 


- 19350/RT 


k f .. - 5.0 x 10 10 T 2 * 80 e - 1400 °/ RT 

I / VX1 


(2.3829$ + ■ - 7 - 6 — 3.5714) 0 574 

k f,viii - ™ ‘ 106 ( ^F" > 


T 3.29 e -16500/RT 


The minimum ranges of validity of the above model are 
0 . 5 < p £ 1 . 0 atm 
0 . 5 < $<4.0 


s 


U 



2% SiH 4 < X < 20% SiH 4 

800 K < T < 1250K 
— O ” 

800K < T < 3000K 

Typical results obtained using the global model are compared 
with results using the compelte mechanism in Table I in figs. 

16-20. Generally, the global model tends to replicate the igni- 
tion phase extremely well. On the other hand, the final approach 
to equilibrium is either too rapid (figs 16 - 18,20) or too slow 
(fig. 19) . However, previous work on global modeling (eg, ref. 5) 
clearly indicates the primary importance of proper replication of 
ignition behavior, in which case the model proposed above should 
yield good-to-excellent results when used in computer codes 
which model complex flows. 

The forthcoming Contractors Report, mentioned earlier, will 
present additional results obtained using the mechanism in Table I, 
including estimates of reaction times, as %/ell as ignition delay 
times . 
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3. EFFECTS OF HYDROGEN ADDITION ON HYDROCARBON IGNITION 
AND FLAME STABILIZATION 


Current interest in the use of hydrocarbon fuels in 
SCRAMJET engines requires that methods be sought for reducing 
the veil-documented lengthy ignition delays and reaction times 
of hydrocarbon fuels. Two possible approaches to reducing 
these times are: (1) injecting (relatively) small quantities 
of hydrogen along with the hydrocarbon fuel with the hone that 
ignition delay times will be more like those for H 2 than for 
the CxHy; and (2) regeneratively heating the hydrocarbon prior 
to its injection into the combustion chamber so as to pyrolyze 
(thermally crack) it, with the expectation that substantial 
quantities of hydrogen will be among the pyrolysis products. 
Studies are underway to investigate both possibilities. Some 
preliminary results are presented here, with additional results 
to be reported upon in a future status report. 

Ignition delay times are shown in fig. 1. Calculations 
have been carried out selecting propane (C^Hg) as the represen- 
tative hydrocarbon fuel using the chemical kinetic mechanism 
for propane in ref. 1. To date, runs have been made using 
H 2 /CgHg fuel mixtures up to o3% (by volume) H 2 . As can be 
seen in fig. 1, at 1.0 atm, $ = 1.0, the reduction in ignition 
delay time compared with that for pure C^Hg is slight. For 
example, at T Q = 1000K, tj^ t 0.05 sec for CgHg and about 
0.025 sec for the 33% H 2 /67% C^Hg mixture. On the other hand, 
t ID = 10 -4 sec for pure H 2 at T Q = 1000K (ref. 2) . It is 
recognized, however, that 33% H 2 in the fuel mixture represents 



very little hydrogen by mass. As an example, the mass equiva- 
lent of the 20% SiH^/80% mixture as applied here would result 
in a fuel mixture of 98.88% Hj/1.12% CjHg by volume. Clearly, 
then, calculations must be carried out for fuel mixtures sub- 
stantially in excess of 33% H 2 (by volume) . A similarly wide 
range of H 2 /CxHy ratios was utilized in the experiments of 
Cookson (ref. 3), shown schematically in fig. 2. His tests, the 
results of which are in fig. 3, would span the range from about 
3% n 2 to nearly 96% H 2 (by volume) had his "main" fuel been 
C^Hg rather than kerosene. It is interesting to note that using 
mode A (see fig. 2), ignition could be achieved at reasonable 
kerosene injection pressure ratios using relatively small quanti- 
ties of H 2 (about 3% to 33% by volume) , whereas mode B required 
much more substantial quantities of H 2 (33% - 96% by volume) . 

Blowout limit correlations obtained using the PSR code 
are shown in figs.. 4 - 7. Fig. 4 duplicates the conditions of 
Longwell and Weiss (ref. 4) and, as can be seen, the calculated 
curve compares well with the experimental one, thereby increasing 
confidence in the chemical kinetic mechanism in ref. 1. In figs. 

5 and 6, for CgHg/H 2 mixtures, little improvement is noted in 
the flame stabilizing properties of these mixtures over pure 
CgHg up to 20% Hj. A similar conclusion holds for CH^/Hj mixtures 
as seen in fig. 7. Additional results will be obtained for 
CxHy/H 2 mixtures having higher concentrations of hydrogen. 

An investigation has been made of the likely equilibrium 
products resulting from propane pyrolysis. The pyrolysis reaction 
was assumed to be (in accordance with ref. 5) 

3o 



9 



§ C 3 H 8 "*■ a C 3 H 8 + b C 2 H 4 + d CH 4 + e H 2 

| ln thlc case ' the equilibrium concentration of H 2 in the ranges 

I * 800K i T i 1200K ' 0.5 < p < 1.0 atm is calculated to be about 

| 40% by volUB ® < see 0)- It is also of interest to note that 

C ' in thcse ran * es ' the remaining 60% by volume is almost entirely 

C 2 H 4 (see fi 9- 9) • Streamtube and PSR calculations with these 
| concentrations as input to assess ignition and reaction times 

■i? flame stabilization behavior are planned. 
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4. INVESTIGATION OF THE RANGES OF APPLICABILITY OF CHEMICAL 
KINETIC MODELS OF HYDROGEN-AIR COMBUSTION 
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This section is an upgraded and refined version of 
Section 4 in our previous Semi-Annual Status Report and 
represents a final report on this subject. 
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4.1. Introduction 


4.1.1 Motivation for the Present Study 

The study of hydrogen-fooled supersonic combustion 
ramjets (scramjets) has been a major part of the research 
programs of the Hypersonic Propulsion Branch at the Langley 
Research Center. Scramjets take advantage of high flignt Macn 
numbers {greater than 5) to achieve propulsion efficiencies 
greater than that of ramjet engines.^ propulsion efficiency 
is a measure of the thrust energy out divided by the comoustion 
energy in. Due to the great flight speeds of scramjets (Mach 
numbers range from 4 to 7), there are time constraints for 
combustion within a combustor of reasonable size. 

As a result of these time constraints# minimal ignition 
delay times are desirable. Figure 1 shows a plot of the log of 
the ignition time versus initial temperature for various fuels. 
Hydrogen fuel has a very low ignition delay time relative to 
hydrocarbon fuels. It is this low ignition delay time which 
makes hydrogen a valuable fuel for scramjets. It becomes a 
necessity, therefore, to study in detail hydrogen-air 
combustion. 

Most studies, thus far, have been experimental 
investigations. This dependence on empirical results is due to 
the complexity of the flow around fuel injectors with 
three-dimensional geometries, which are not easily treated 
analytically. Numerical solutions have generally been 
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restricted to two- or three-dimensional parabolic flow with 
oversimplified chemistry models of the H2~air system. These 
numerical solution schemes are applicable only in the parabolic 
flow region veil downstream of the disturbance 



To 


FIGURE 1 IGNITION t)ELAY TII-.E VS. 

INITIAL TEI-.FERa.TURE 

caused by the transverse fuel injection used by scramget 
combustors in order to achieve rapid mixing and reaction. A 
priori knowledge of the extent of fuel mixing, ignition, and 
reaction is required to initiate calculations. 

In spite of the difficulty in obtaining quantitative 
information, a sufficient data base has been established to 
define a scramjet engine concept and to permit fabrication of 
subscale engine models with integrated inlet., combustor, and 
nozzle components. The current scramjet is designed to operate 

- 4f- 







at stagnation temperatures between 900*K and 2200 *r 2 
which correspond to a flight Mach number range from 4 to 
7. In ground tests of subscale engine models^, however, 
problems were encountered in obtaining ignition and 
sustaining reaction at test conditions where ignition ana 
sustained reaction were expected. 

The need to better understand the chemical mechanism 
of the ignition and reaction of H 2 ~air mixtures at 
conditions typical of a scramjet combustor has lea to many 
analytical studies. Most of the analytical studies at the 
Langley Research Center, as well as this present stuay, 
use a computer program (references 4 and 5) to solve 
flowing, chemical-kinetic, isobaric, stream-tube problems 
involving many chemical species. It is known that the 
computational time requirements for any computer program 
employing detailed chemical kinetics is proportional to 
the number of species and reactions being treated. It is 
the purpose of this effort to develop a method which 
reduces these numbers in the study of H 2 - air combustion 
and at the same time preserve the correct physical- 
chemical behavior. The result of this work will be to 
reduce computer time and computer storage requirements. 
Methods for reducing computational times are presented in 
section two of this work. 
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4. 2. Methods to Reduce Computational Times 


The method proposed in this work to reduce 
computer run times is an extension and refinement of 
the method proposed by Chirrtz in reference 9. This 
method involves the tracking of a "trigger" species 
which is used to determine whether the flow is in an 
"ignition" mode or in a "combustion" mode. Ignition 
and combustion are defined in the classical sense. 
Ignition delay time is taken to be the time required 
for the temperature increase to reach five percent of 
the overall temperature increase: 

T ig • - T 0 + 0. 05 (T e q - T 0 ) ( 1 ) 

When the temperature is less than the ignition 
temperature, the flow is in the "ignition” mode. 
Otherwise, the flow is in the "combustion" mode. 

(See Figure 2). 

It will be shown that while an extensive 
chemical package is needed to describe the "ignition" 
mode, a smaller package is sufficient to deal with 
the "combustion" mode. Evans and Schexnayder 1 ^ 
concluded that a 25-reaction scheme involving 12 
species (designated 25(12) herein), was required to 
describe "ignition" processes, while an 8(7) sufficed 
to deal with the "combustion" mode. The work 
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discussed herein describes efforts to minimize the 
computational time requirements once the "combustion" 
process is initiated (i.e. a smaller package or one 
that takes less computational time than the b(7) 
system might suffice to deal with the "combustion" 
mode) . 

A 37(13) system is used as the "test" 
chemical-kinetic package (Table 1). Any smaller 
system is "tested" against the 37(13) system to see 
if the correct physical chemical behavior is 
preserved. This is accomplished in the following 
way* 

First, the full 37(13) system is run from tne 
initial state to equilibrium. Next, the test system 
is run from the point where the "combustion" in the 
37(13) system initiates. The test system's initial 
values (pressure, temperature, equivalence ratio ana 
chemical composition) are that of the 37(13) system's 
values at the point of "combustion". Temperature- 
time profiles and chemical behavior are compared (see 
Figure 3). If the difference between the 37(13) 
system and the test system is within an acceptable 
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range, the test system will be sufficient to deal 
with the "combustion" process. 

Graphs of the mass fraction of the trigger 
species versus time and ignition temperature versus 
mass fraction of the trigger species at ignition are 
made in accordance with reference 9. The cases from 
the 37(13) system serve as the data base. 

An 8(7) system (Table 2) is tested to confirm 
the results of Evans and Schexnayder . The 
8(7) system is run from the initial state to 
equilibrium to determine if there are any initial 
conditions where the 8(7) system describes the entire 
process (i.e. both ignition and combustion). 

In order to further reduce computational 
times, a 2(5) global model (Table 3) developed oy 
Rogers and Chinitz in reference 11 is tested. The 
method proposed in that paper must be altered to 
include the effects of pressure. The Arrhenius 
equation will have the form: 

kfi « Aj( , P )T Ni exp(-Ei/RT) (2 ) 

The values of the parameters may be different from 
those in reference 11 because they are fixed 
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arbitrarily so that the 2(5) system describes the 
37(13) system. 

Lastly# the so-called "partial equilibrium" 
assumption of reference 12 is examined in connection 
with an 8(7) system in an effort to minimize 
computational times. The 8(7) mechanism consists of 
the following reactions: 


H 2 0 + 

0 

*e 

OH 

+ OH 

(i) 

H 2 0 + 

H 

- 

OH 

+ h 2 

(ii ) 

°2 + 

H 

a 

OH 

+ 0 

( iii ) 

h 2 + 

0 

« 

OH 

+ H 

(xv) 

h 2 0 + 

M 

ct 

H 

+ OH 

(v) 

H 2 ♦ 

M 

a 

H 

+ H + M 

(vi) 

0 2 + M 

z= 

0 

+ 0 + M 

(vii ) 

OH + 

M 

B 

0 

+ H + M 

(vxii) 


During the ignition delay period, all eight 
reactions play a role; however, in the "combustion" 
mode, the bimolccular shuffle reactions occur go 
rapidly in both directions that under some conditions 
they may be basically in equilibrium. The "partial 
equilibrium" assumption takes ractions (i)-(iv) to De 


in equilibrium (infinitely fast reaction rates), 
while reactions (v)-(viii) continue to occur at a 
finite rate. 

Infinitely fast reaction rates can be app- 
roximated by assigning very large values to the 
preexponential factor in the forward reaction rate 
constant (e.g. 10^®). The backward **eaction rate 
constant is then very large to satisfy k = k f /k b 
The values of the preexponential factors of reaction 
(i) to (iv) are determined so that the system des- 
cribes the "combustion" process. 

To obtain these results, 45 one-dimensional, 
constant pressure, ^“air computations were 
performed .in the ranges: 

0.5 < p < 1.5 

850 < T <1200 

0.5 < p <1.0 

where p is the equivalence ratio, T is the 

temperature in *K, *nd p is the pressure in atmos- 
pheres. A listing cf these individual cases is given 
in Table 4. Comparison of the computational times 
are then made and recommendations are discussed. 
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4.3. Results end Discussion 



4.3.1 Summary 

All 45 cases in Table 4, as well as the nine 
cases narked with an exclamation point, were 
performed on the 37(13) and 8(7) systems, 
respectively, in both the ignition and combustion 
modes. In addition, the 15 cases marked with an 
asterisk were tested on the 8(7) system, the gloDal 
model, and the "partial equilibrium" model starting 
at the ignition point. These 15 cases were selected 
to encompass the entire range of interest. 

What follows are the results of these test cases. 
First, the results from the 37(13) system including 
the "trigger" concept are presented. Next, the 8(7) 
system in its entirety and in the combustion mode 
alone is discussed. Global and partial equilibrium 
results in the combustion mode follow. Lastly, 
temperature-time profiles and computational times are 
compared for all four systems. 

4.3.2 The 37(13) System 

4*3.2. 1 Temperature-Time Profiles 

45 cases were run to serve as a basis 
of comparison for the test systems described 
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previously. Sufficient data points were obtained for 
all cases and each case yielded a characteristic 
S-shaped curve for the temperature versus time 
profile. Por each case, it was necessary to obtain a 
data point at the ignition temperature, defined in 
the classical sense as 5% of the temperature rise 
from the initial state to the final equilibrium 
state, in order to set initial conditions for the 
test cases. Problems arose, however, when trying to 
pinpoint the ignition temperature. As a result, 
satisfactory ignition conditions could not be 
calculated for nine cases due to time constraints. 
These cases were omitted from data base and the 
remaining 36 cases are given in Table 5 with 
corresponding ignition temperature, time, and mass 
fraction of OH. 

Tables 4 and 5 show that the ignition delay time 
is a function of initial temperature, pressure, and 
equivalence ratio. This function is complicated 
since it is inversely proportional to initial 
temperature, nearly independent of equivalence ratio 
except at low temperatures (850*K and 900*K), 
proportional to pressure at temperatures of 850*K ana 
900*K, and is inversely proportional to pressure at 
temperatures greater than or equal to 10Q0*K in the 
ranges studied. These recults are similar to the 
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results reported by Rogers and Schexnayder (reference 
2) in their extensive ctudy which included 6u 
reactions and 20 species. Since the results are 
similar to those of Rogers and Schexnayder, the 
37(13) system serves as a good basis for the 
remaining systems to be compared with. 

4. 3. 2. 2 The "Trigger" Concept 

In order to test systems in the 
"combustion" mode it is necessary to determine when 
the transition to the "combustion" mode takes place. 
It was to this end that the concept of a "trigger" 
species was first proposed by Chinitz in reference 9. 
The concept involes the tracking of the mass fraction 
of one of the species. When its value becomes 
greater than a preset value, the system is said t.c be 
in the "combustion" mode. 

There are several requirements for a species to 
serve as an effective trigger from ignition to 
combustion. The concentration of an effective 
trigger species must be single-valued (the value cf 
the mass fraction at the ignition point is not again 
obtained) and undergo a sufficiently large change 
during the ignition process that its crossover into 
the combustion mode is unmistakenble. It was 
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It is of interest to estimate the trigger point 
for all values of equivalence ratio, pressure, and 
temperature studied. A plot of the ignition 
temperature vs the mass fraction for OH at the 
ignition point is shown in Figure 5. A least-squares 
linear f . : of In Tjg is Yo H ^g was determined 
to be 

In T ig - 0.085 In Y 0Hig + 7.65 ( 3 ) 

This linear relationship predicts values within a 
factor of 2.5 (generally well within) for all 
pressures an' equivalence ratios examined. This 
provides a large advantage over the non-linear 
relat ionsh ip. developed in reference 9 which was 
limited to an equivalence ratio of 1. Kith the 
selection of a triyger species whose value can be 
estimated, the test systems can now be examined. 

— y ^ — 




4.3.3 The 8(7) System 
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4. 3. 3.1 Ignition and Combustion 

The 8(7) sy totem was run from time 
zero to equilibrium for the nine cases marked with an 
exclamation point in Table 4. Five of these cases 
(3* 6* 9, 12* 15) represent a set of constant 
pressure* constant equivalence ratio conditions with 
varying initial temperature. These were selected to 
determine how the initial temperature affects the 
ignition delay time when compared to the 37(13) 
system. The remaining four cases (13, 22, 25, 29 ) 
were arbitrarily picked in an attempt to see the 
effects of pressure and equivalence ratio. 

Figures 6-10 show a comparison of temperature- 
time histories for the 37(13) and 8(7) systems for 
cases 3, 6, 9, 12, 15, respectively, in both the 
ignition and combustion modes. Figures 6 and 7 show 
the delay time to be more than an order of magnitude 
lower for the 8(7) system at low temperature. 

Figures 8-10, however, show that the delay time for 
the 8(7) system is greater than the 37(13) system at 
higher temperatures. It is therefore concluded that 
a transition takes place between an initial 
temperature of 900*K and 1000*K. (It should be noted 
that a ten percent mixture of hydrogen in air has an 

i 3 

ignition temperature of between 893*K and 1020*K ). 
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Pigurc 11 shove a plot of ignition delay time vs. 
1000T o “' for the 8(7) and 37(13) systems at a 
pressure of one atmosphere and an equivalence ratio 
of 0.5. The two curves cross at a temperature of 
970*K where their ignition delay tines are equal. 

These results imply tnat for any constant 
pressure and constant $ there is a T 0 for which 
the two systems have the same ignition delay time. 

At these conditions, the 8(7) system can descnce the 
37(13) system in both the ignition and combustion 
modes. It also might be expected that at 
temperatures away from this particular initial 
temperature, the 8(7) system might lag behind or 
proceed ahead of the 37(13) system when testea at 
ignition conditions in the "combustion" mode. 

The remaining four cases all had initial 
temperatures at or above 1000*K. The results of 
these cases were again that the delay time was higher 
for the 8(7) system than that of the 37(13) system. 

No significant effects of pressure or equivalence 
ratio were found in these limited case. The 8(7) 
system was now tested in the "combustion" mode witn a 
focus on the effect of initial temperature. 



4. 3.3.2 The "Combustion" Mode 


The 15 esses narked with on asterisk 
in Table 4 were run for the 0(7) system in the 
"combustion node." Por all pressures, temperatures, 
and equivalence ratios studied, the 8(7) system haa a 
higher equilibrium temperature than that of the 
37(13) system. This result was expected in that the 
8(7) system neglects nitrogen dissociation. It was 
also found that for all conditions, the 8(7) system 
lagged behind the 37(13) system after the ignition 
point then equaled and surpassed the 37(13) curve 
(see Figure 12). This result was contrary to the 
expectations of section 4. 3. 2.1. It was thought tnat 
for some initial temperatures the 8(7) system would 
reach equilibrium more equickly than the 37(13) 
system. This, however, was not the case. 

The amount of lag, however, was found to be a 
function of the relationship between the actual 
ignition temperature and the ignition temperature 
calculated from equation 1 (designated herein as 
"classical" ignition temperature). If the ignition 
temperature used (and also all other initial 
conditions) was between the initial temperature and 
the classical ignition temperature, then the 8(7) 
system described the 37(13) system very well. If the 
temperature used was greater than the classical 
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ignition temperature, there was an appreciable lag 
(see Figure 12). Therefore, as the difference 
between the initial temperature and ignition 
temperature used increased, the lag became greater. 

No cases were found where the 0(7) system reached 
} equilibrium more quickly than the 37(13) system. It 

i 

is concluded, therefore, that ignition temperatures 
lower than the classical ignition temperature should 
be used. It was not determined what percent below 

i 

five percent should be used. 

4,3.4 The Global Model 

The global model in Table 3 was first 
proposed in reference 11. The reaction rate constant 
for each reaction is adjusted so that the global 
model describes the 37(13) system in the comoustion 
mode. 

Before determining these adjusted rate 
constants, it was useful to observe how varying each 
rate constant affected the temperature-time profile. 
It was observed that increasing (the 

subscript 4 refers to the first reaction; subscript 5 
refers to the second reaction using the notation in 
reference 11) slows down the reaction. This was 
referred to as a "lagging" system ir. the previous 
section. Increasing k^ speeds up the reaction. 

Tne syotera is very sensitive to an increase or 
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decrease in k^g while being rather insensitive to 
a variation in k £4 * It was decided to fix Ni and 
Ei at the came values as in reference 10; namely , 

E 4 * 4865 cal/mol 
N 4 0 -10 

E 5 “ 42,500 cal/mol 
N 5 ■ -13 

Due to the insensitivity of the system to 
k £ 4 , it was also decided that A 4 could be 
retained as in reference 11 ; 

A 4 U) = (8.917**. + 31.433)<f> - 28.950) x 10 47 , cm 3 /mol-s. 


It should be noted that A 4 is a function of jS only. 

The effect of pressure, which was neglected in 
reference 11 , was included in the determination of 
Aj. As a result, A 5 is a function of pressure 
and equivalence ratio and has the form A s (^,p). 

The sensitivity of the system to kf 5 permitted 
substantial changes to be made to the temperature 
profiles produced by the global system. 

Table 8 shows the values of A 4 and A 5 
determined in order to describe the temperature-time 
profile required. Figure 13 shows a plot of A 4 
(which is independent of pressure) vs. ^ . Figure 14 

shows plots of A 5 vs. pressure for various (j) 's. 
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As other authors have reported , large discrepancies 
exist between results for equivalence ratios less 
than or greater than one. The relationship between 
A5 and pressure and equivalence ratio must be 
broken up into two equations for the ranges studied: 

F ° r ♦ 1 1 

ft 5 .(8.80 + 5.85/$ - 3.67$ - 4.80p - 3.20p/$.+ 2.00p-$j MO 
For d < 1 

A . O^p' 1 V°- 7 (8.80 + 5.85/$ - 4.80p - 3.20p/$ ♦ 2.00p$)« 
® . 10 

These equations predict values of A 5 to within 
10 % for the ^ 's studied.. 

In all cases, the global mode lags behind the 
37(13) system for a short period following ignition. 
Then the global model equals and exceeds the 37(13) 
system and reaches equilibrium (see Figure 15). The 
equilibrium temperature of the global model should be 
higher than that of the 37(13) system due to 
dissociation . 

This again, however, is controlled by the 
ignition temperature used by the global model. If 
the ignition temperature used is greater than the 
classical ignition temperature, then the equilibrium 
temperature is less than the classical ignition 
temperature. If the ignition temperature is less 



than the classical , then the equilibrium temperature 
reached is greater than the 37(13) system as 
expected. 

Initial temperature also has an effect on the 
accuracy of the global model. As the initial 
temperature increases, the difference between the 
global model and the 37(13) system increases, ho 
effect on pressure or equivalence ratio was 
found. 

The global model predicts very well the 
temperature-time profile of the 37(13) system. These 
results, however, show the importance of the ignition 
temperature used. It is suggested that ignition 
temperature be a variable to be tested in future 
work. Nonetheless, the global model shows promise in 
accurately reproducing temperature- time profiles with 
a minimum of computational time and computer storage 
requirements . 

4.3.5 The "Partial Equilibrium" Assumption 

The "partial equilibrium 3 assumption is 
already in use in combustion analysis; however its 
applicability to this present Gtudy in the ranges of 
interest has not been determined. This assumption 
states that the bimolecular shuffling reactions occur 


bo rapidly that the ere basically in equilibrium. If 
this is the case, the kinetics equations can be 
replaced by the algebraic laws of mass action as 
discussed in reference 12. Therefore, a number of 
partial differential equations are replaced by 
algebraic equations and the remaining partial 
differential equations are simplified. If 
appropriate criteria are specified as to when this 
assumption is applicable, computer running times will 
be greatly reduced. 

It was to this end that an approximation to tne 
"partial equilibrium" assumption was tried. Instead 
of incorporating algebraic equations into the 
computer program, infinitely fast reaction rates 
would be approximated by using extremely large 
numbers (on the order of 10^) for the 
preexponential factor. It was anticipated that this 
would sufficiently approx ; mate the "partial 
equilibrium" assumption. 

Substantial numerical difficulties arose in 
trying the approximation. First, the largest 
preexponential factors that could be used were of the 
order of 10^. Still, this was at least 5 orders 
of magnitude greater than any other preexponential 
factors, Hext, the precision of the program needed to 
be upgraded two orders of magnitude by adjusting tne 



EMAX parameter. Pinally, this approach produced 
erroneous results ouch as temperatures exceeding the 
equilibrium temperature when the Gear integration 
procedure was used. The slower Adams method was then 
used which gave results which were physically 
plausible. Due to these problems, running times were 
very long. This was not the main concern, however. 
Rather, the primary goal was to test the partial 
equilibrium assumption against the 37(13) system in 
the ranges studied. 

Figures 16-30 show how the approximation to the 
"partial equilibrium" assumption compared to the 
37(13) system. For low temperatures (850*K, 900*K) 
and for all cases with an equivalence ratio of 0.5, 
the partial equilibrium assumption reproduced the 
37(13) curve well. In the remaining cases, however, 
the assumption did not approximate the profile 
accurately. The shape of the curve was not even 
preserved . 

It appears that the approach to approximating 
the "partial equilibrium" assumption has a narrow 
region of validity; namely, low initial temperatures 
and low equivalence ratios. As stated before, no 
conclusion can be drawn about the assumption itself. 
It is suggested that future work include the 
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incorporation of algebraic laws of mass action into 
the computer program to nore precisely determine the 
validity of the assumption. 

4.3.6 Comparison of the Pour Systems 

4,3. 6.1 Temperature-Time Profiles 

Figures 16-30 show a comparison 
of the 37(13) systems for the 15 cases marked with an 
asterisk in Table 4. The ignition mode of all 15 
cases is that of the 37(13) system which continues in 
the combustion mode. The test systems start at the 
ignition point and reach an equilibrium temperature 
different from the 37(13) system. 

The lag of the 8(7) system is shown in these 
figures. The lag is more a function of the 
relationship between the ignition temperature used 
and the classical ignition temperature than the 
initial temperature, pressure, or equivalence ratio. 
The 8(7) system reproduces the temperature-time curve 
very well, as expected. 

The inital lag and eventual higher equilibrium 
temperature of the global model is shown in all 15 
cases. Initial temperature has a substantial effect 
on the accuracy of the global mode. As the initial 
temperature increases, the temperature differences 
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between the global code! and the 37(13) carves become 
large. These perhaps could be corrected by starting 
the global model at a different ignition temperature. 

Partial equilibrium curves accurately reproduced 
the temperature-time curves for cases 2, 4, 9 , 11, 

15, 18, 20 and 31. This suggests that the partial 
equilibrium assumption might suffice at lower initial 
temperatures and equivalence ratios. The numerical 
difficulty in obtaining these results sheds 
uncertainty on this conclusion. 

Figures 31 and 32 represent the results of using 
an ignition temperature higher than the classical one 
to cases 11 and 15, respectively. The 8(7) system is 
observed to lag behind the 37(13) system to a greater 
extent than usual. The global model, in fact, has an 
equilibrium temperature which is 200 degrees below 
normal. These figures reinforce the importance in 
the selection of initial conditions for any test 
system. 

4. 3.6.2 Computational Time 

Results of the average 

computational times for the various systems are shown 
in Table 9. The global model has the lowest 
computational time as expected, with an average time 






of epproximatley 1/0 that of the 37(13) Bystem. 
Excessive computational time of the partial 
equilibrium system was related to the relatively nigh 
reaction rate constants employed, the use of the 
Adams method rather than the Gear method, and the 
precision which was required to perform "partial 
equilibrium" calculations. 




4,4. Conclusions and Recommendations 
4.4.1 Summary 

A data base* in the ranges of initial 
temperature between 850*K and 1200*K, pressure 
between 0.5 and 1.00 atmospheres and equivalence 
ratio between 0.5 and 1.5, was established with the 
37(13) system serving as the basis. Due to the large 
amount of data produced, it is necessary to clarify 
these results. First, these results are summarized. 
Next, conclusions are enumerated. Finally, 
recommendations for future work are given. 

The first part of the work was selecting a good 
trigger species. It was concluded that the hydroxyl 
radical, OH, serves as an effective trigger. The 
value of the mass fraction of OH at ignition can De 
estimated in the ranges studied within a factor of 
2.5 (generally well within) according to equation 
32. 

Next the 8(7) system was studied in the 
"ignition" and "combustion" modes together and then 
in the "combustion" mode alone. It was concluded 
that a transition takes place between an initial 
temperature of 900*K and 1000*K (at a pressure of one 
atmosphere and an equivalence ratio of 0.5) where the 
ignition delay time of the 8(7) system switches from 
lesG than to 
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greater than the 37(13) system. It was further 
concluded that at an equivalence ratio of 0.5* 
pressure of 1.0 atmosphere and initial temperature of 
970+30*K the 8(7) system can approximate the 
temperature-time profiles of the 37(13} system in 
both the ignition and combustion modes. It is 
recommended that future work include the 
investigation of conditions where the 8(7) system can 
describe the 37(13) system in both inodes. 

In the combustion mode alone, the 8(7) system 
reproduces the temperature-time profiles very well 
for all cases studied. The variables which were 
major factors in the accuracy of these profiles were 
the temperatures and species' concentrations used at 
ignition. If the temperature used was less than, but 
close to, the classical ignition temperature, then 
the profiles v.,re very accurate. If, however, the 
temperature used was greater *han the classical 
ignition temperatures, then the 6(7) system lagged 
behind the 37(13) system somewhat. Future work, 
should include the investigation of this phenomenon. 

Next, the global model was studied in the 
combustion mode. Preexponential factors were 
determined with a rather large discrepancy between 
results for an equivalence ratio less than or greater 
than one. For all cases, the global model lags 
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behind the 37(13) eyotera initially then exceeds it at 
equilibrium Two variables had an effect on the 
accuracy of this codel. First, the conditions used 
at ignition had a substantial effect on accuracy 
again* A lag was produced by using a temperature 
greater than the classical ignition temperature. 
Secondly, as the initial temperature increases the 
accuracy of the global model decreases. 

Finally, an approximation to the "partial 
equilibrium" assumption was tried. Numerical 
difficulties arose in the execution of this approxi- 
mation. Nonetheless, it was concluded that at low 
temperatures and low equivalence ratios the approach 
to approximating the partial equilibrium assumption 
is valid. No conclusions could be drawn about the 
"partial equilibrium" assumption itself. 

Results for the various systems show 
the global model to have the lowest average 
computational time followed by the 8(7) system, the 
37(13) system, and the “partial equilibrium" model. 
The two-step global model has a computational time of 
approximatley 1/8 that of the 37(13) system. 

4.4.2 Conclusions 

The following conclusions can be drawn based on 
the results of this study: 

- C.? - 
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(1) The hydrosyl radical, CH, serves as an 
effective trigger species. 

(2) A transition takes place between 900*K and 
1000*K (at a pressure of 1 atmosphere and 
an equivalence ratio of 0.5) where the 
ignition delay tine of the 8(7) system 
switches from less than to greater than the 
37(13) system. 

(3) At a pressure of one atmosphere, an 
equivalence ratio of 0.5 and an initial 
temperature of 970+30*K the 8(7) system can 
accurately approximate the temperature- 
time profile of the 37(13) system in both 
the ignition and combustion modes. 

(4) For all cases, the 8(7) system describes 
the 37(13) system very well in the 
combustion mode. 

(5) A lag is produced in the 8(7) and global 
models when using ignition temperatures 
greater than the classical ignition 
temperature . 

(6) The global model satisfactorily reproduces 
the temperature-time profiles of the 37(13) 
system with a minimum of computational time 
and computer storage requirements. 

(7) As the initial temperature increases, the 
accuracy of the global model decreases. 

(8) No conclusions can be drawn concerning the 
"partial equilibrium" assumption. 

(9) The approach to approximating the "partial 
equilibrium" assumption has a narrow region 
of validity? namely, low initial 
temperatures and low equivalence ratios. 


7 ° 



(10) The global model can approximate the 37(13) 
system with an average of 1/0 the 
computational tine. 

4.3 Recommendations 

Based on the results of this study, the 
following recommendations are nade for future worx: 

(1) Investigate conditions where the 0(7) 
system can accurately describe both the 
ignition and combustion processes. 

(2) Study the effects of varying ignition 
conditions; namely, above and below the 
classicial ignition temperature. This will 
show the relationship between ignition 
conditions and the lag of the system. 

(3) Continue to investigate the two-step glooal 
model as it seems to have the most promise 
in reproducing temperature-time profiles 
with a minimum of computational time. 

(4) Incorporate the algebraic laws of mass 
action into the computer program to more 
precisely determine the validity of the 
"partial equilibrium" assumption. 
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Table 1 - 37(13) System 


REACTION 

NUMBER 


1 

M 

4 

1*02 

O 

1*0 

4 

1*0 

2 

M 

4 

1*H2 

m 

1*H 

4 

1 * H 

3* 

M 

4 

1*U20 

m 

1*H 

4 

l*0n 

4 

1*H 

4 

1*02 

m 

1*H02 

4 

M 

5 

N 

4 

1*N02 

m 

1*N0 

4 

1*0 

6 

M 

4 

1*N0 

SR 

1*N 

4 

1*0 

7 

H 

4 

1*H202 

B 

1*0H 

4 

l*0rt 

8 

M 

4 

1*03 

e 

1*02 

4 

1*0 

9 

1*0 

4 

1*H 

B 

1*0H 

4 

M 

10 

1*H20 

4* 

1*0 

m 

1*0H 

4 

l*0ri 

11 

1*H2 

4 

1*0H 

m 

1*H20 

4 

1 *H 

12 

1*02 

4 

1*H 

B 

1*0H 

4 

1*0 

13 

1*H2 

4 

1*0 

K 

1*0H 

4 

1 *H 

14 

1*H2 

4 

1*02 

MX 

1*0H 

4 

l*OU 

IS 

1*H 

4 

1*H02 

B 

1*H2 

4 

1*02 

16 

1*H2 

4* 

1*02 

B 

1*H20 

4 

1*0 

17 

1*H 

4 

1*H02 

a 

1*0H 

4 

l*Qri 

18 

1*H20 

4 

1*0 

B 

1*H 

4 

1*H0/ 

19 

1*0 

4 

1*H02 

MS 

1*0H 

4 

1*02 

20 

l*OH 

4 

1*H02 

B 

1*02 

4 

1*M20 

21 

1*H2 

4 

1*H02 

B 

1*H20 

4 

l*0tl 

22 

1*H02 

4 

1*H2 

B 

1*H 

4 

l*h202 

23 

1*H202 

4 

1*H 

R 

1*0H 

4 

1*«120 

24 

1*H02 

4 

1*0H 

C 

1*0 

4 

1*11202 

25 

1*H02 

* 4 

1*H20 

s 

1*0H 

4 

l*n202 

26 

1*H02 

4 

1*H02 

m 

1*H202 

4 

1*02 

27 

1*0 

4 

1*03 

= 

1*02 

4 

1*02 

28 

1*03 

4 

1*H0 

c 

1*N02 

4 

1*02 

29 

1*03 

4 

1*H 

B 

1*0H 

4 

1*02 

30 

1*03 

4 

1*0H 

MS 

1*02 

4 

1*ho2 

31 

1*0 

4 

i*N2 

B 

1*N0 

4 

1 * N 

32 

1*H 

4 

1*N0 

m 

1*0H 

4 

1*N 

33 

1*0 

4 

1*N0 

c 

1*02 

4 

1*N 

34 

1*N02 

4 

1 *H 

ss 

1*N0 

4 

1 *0H 

35 

1*N02 

4 

1*0 

= 

1 *N0 

4 

1*02 

36 

1*H02 

4 

1*N0 

B 

1*20 

4 

l*0n 

37 

1*03 

4 

1*H02 

MS 

2*02 

4 

l*0tl 


^2 is an inert species 
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TABLE 2 


8(7) SYSTEM 


REACTION 

NUMBER 




REACTION 



1 

H 


02 

at 

0 

4* 

O+K 

2 

M 

+ 

H2 

e 

H 

4 

H+K 

3 

M 


H20 

at 

H 

4 

OH-v 

4 

M+0 

+ 

H 

•t 

OH 

4 

K 

5 

H20 

4 

0 

* 

OH 

4 

H 

6 

H2 

4 

OH 

v» 

H20 

4 

H 

7 

02 

4 

H 

it 

OH 

4 

0 

8 

H2 

4 

0 

e 

OH 

4 

H 


N'2 is an inert species 



TABLE 3 


STEP GLOBAL MODEL 


- ? 




H2 + 02 - 2*0H 

(1) 

k fc 

2*0H + H2 -2*H20 

(2) 

N_ is an inert species 
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CASES 

TABLE 4 
$ 

- 45 CASES 

p, atm 

T, K 

1 

0.5 

0.5 

850 

* 2 

0.5 

0.75 

850 

! 3 

0.5 

1.0 

850 

* 4 

0.5 

0.5 

900 

5 

0.5 

0.75 

900 

! 6 

0.5 

1.0 

900 

7 

0.5 

0.5 

1000 

8 

0.5 

0.75 

1000 

'!* 9 

0.5 

1.0 

1000 

10 

0.5 

0.5 

1100 

*11 

0.5 

0.75 

1100 

112 

0.5 

1.0 

1100 

113 

0.5 

0.5 

1200 

14 

0.5 

0.75 

1200 

*115 

0.5 

1.0 

1200 

16 

1.0 

0.5 

850 

17 

1.0 

0.75 

850 

*18 

1.0 

1.0 

850 

19 

1.0 

0.5 

900 

*20 

1.0 

0.75 

900 

21 

1.0 

1.0 

900 

122 

1.0 

0.5 

1000 

*23 

1.0 

0.75 

1000 

• 

ro 

1.0 

1.0 

1000 

*125 

1.0 

0.5 

.1100 

26 

1.0 

0.75 

1100 

27 

1.0 

1.0 

1100 

*28 

1.0 

0.5 

1200 

129 

1.0 

0.75 

1200 

30 

1.0 

1.0 

1200 

*31 

1.5 

0.5 

850 

32 

1.5 

0.75 

850 

33 

1.5 

1.0 

850 

34 

1.5 

0.5 

900 

35 

1.5 

0.75 

900 

*36 

1.5 

1.0 

900 

37 

1.5 

0.5 

1000 
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TABLE 4 (Cont.) 


CASES 

♦ 

p, atm 

T, K 

38 

1.5 

0.75 

1000 

*39 

1.5 

1.0 

1000 

40 

1.5 

0.5 

1100 

*41 

1.5 

0.75 

1100 

42 

1.5 

1.0 

1100 

*43 

1.5 

0.5 

1200 

44 

1.5 

0.75 

1200 

45 

1.5 

1.0 

1200 
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TABLE 5 


37(13) DATA 


2 


918 

4 


961 

• 5 

• 


960 

6 


974 

7 


1045 

8 


1053 

9 


1024 

-11 


1160 

13 


1260 

15 


1260 

18 


935 

19 


970 

20 


985 

21 


966 

22 


1027 

23 


1012 

24 


1039 

25 


1170 

26 


1117 

27 


1170 

28 


1275 

29 


1275 

30. 

* - 

1270 

31 


926 

32 


892 

34 


934 

35 


968 

36 


954 

37 


1012 

. 33- 


1024 

39 


1052 

40 


1160 

41 


1163 

42 


1160 

43 


1275 

45 


1303 


lia 


0.48E0 

0.66E-4 

0.80E-3 

0.10E-3 

0.13E-2 

0.13E-3 

0. 1 IE- 1 

0.16E-3 

0.21E-3 

0.51E-3 

0.15E-3 

0.60E-3 

0.12E-3 

0.15E-3 

0.69E-4 

0.90E-3 

0.60E-4 

0.34E-2 

0.29E-4 

0.22E-2 

0.46E0 

0.19E-3 

0.76E-3 

0.40E-3 

0.12E-2 

0.19E-3 

0.48E-2 

0.10E-3 

0.20E-3 

0.20E-3 

0.14E-2 

0.11E-3 

0.12E-3 

0.19E-3 

0.10E-3 

0.16E-3 

0.60E-4 

0.27E-3 

0.45E-4 

0.90E-3 

0.55E-4 

0.30E-2 

0.38E-4 

0.30E-2 

0.35E-4 

0.30E-2 

0. 1 IE- 1 

0.80E-4 

0.27E0 

0.17E-4 

0.80E-3 

0.57E-4 

0.11E-2 

0.11E-3 

0.35E-2 

0.61E-4 

0.20E-3 

0.67E-4 

0.15E-3 

0.11E-3 

0.13E-3 

0.26E-3 

0.10E-3 

0.91E-3 

0.70E-4 

0.91E-3 

0.45E-4 

0.90E-3 

0.59E-4 

0.17E-2 

0.30E-4 

0.29E-2 
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TABLE 6 


- TRIGGER SPECIES vs. TIME 



TIME 

0 

0.1E-5 
0.2E-4 
0.3E-4 
0.4E-4 
0.6E-4 
0.7E-4 
0.8E-4 
0.9E-4 
CUE- 3 
0.12E-3 
*0.14E-3 
0.2E-3 
0.3E-3 
0.4E-3 
0.5E-3 
0.8E-3 


0 

n.68E-9 

0.38E-8 

0.81E-8 

0.17E-7 

0.75E-7 

0.16E-6 

0.35E-6 

0.82E-6 

0.21E-5 

0.15E-4 

n.UE-3 

0.15E-1 

0.19E-1 

0.19E-1 

0.18E-1 

0.17E-1 


ignition point 



TABLE 7 

COMPARISON OF IGNITION TIMES FOR 37(13) AND 8(7) SYSTEMS 


la 

37(13) t ig (sec) 

8(7) Tiq UtO 

850 

6.7 x 10* 1 

1.1 x 10‘ 3 

900 

1.1 x 10“ 2 

6.6 x 10' 4 

1000 

1.3 x 10‘ 4 

2.6 x 10" 4 

1100 

0.5 x 10' 4 

1.5 x 10* 4 

1200 

0.3 x 10" 4 

0.7 x 10” 4 





TABLE 8 

PREEXPONENTIAL PACTORS POR GLOBAL MODEL 


CASE 

A 4 

A 5 

2 

3.8E48 

4.1E65 

4 

3.8E48 

4.4E65 

9 

, 3.8E48 

2.3E65 

11 

3.8E48 

2.3E65 

15 

3.8E48 

1.1E65 

18 

1.2E48 

5.5E64 

20 

1.2E48 

9.0E64 

23 

1.2E48 

6.0E64 

25 

1.2E48 

9.0E64 

28 

1.2E48 

8.0E64 

31 

5.0E47 

6.3E64 

36 

5.0E47 

3.7E64 

39 

5.0E47 

4.0E64 

41 

5.0E47 

4.6E64 

43 

5.0E47 

5.1E64 


CO 



TABLE 9 

COMPARISON OP COMPUTATIONAL TIMES 


SYSTEM 

AVERAGE 

COMPUTIONAL 

TIME (seconds) 

37(13) 

87.5 

8.(7) 

12.3 

Global 

11.4 

Partial 

1258 
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5. APPLICATION OF THE CHARNAL CODE TO TURBULENT REACTION 
RATE MODELING 

Work has been continuing dealing with the use of the 
CHARNAL computer code to assess and refine the turbulent reaction 
rate model discussed in previous status reports. At this time, 
three experiments are being examined using the code: (1) the 

Beach experiment (initial results are contained in a forthcoming 
NASA TM) ; (2) the Kent and Bilger experiment (see fig. 1); 

(3) the Northam experiment (see fig. 1). Some preliminary results 
obtained using the code in its non-finite-rate reaction modes 
are shown in figs. 2-5. In fig. 2, the experimental axial center- 
line temperature profile is shown along with CHARNAL results run 
as free jet with specified axial pressure gradient. As can be 
seen, the assumption of complete reaction results in too rapid a 
temperature increase and too high a maximum temperature. The 
eddy breakup assumption produces opposite results, indicating too 
slow a reaction rate. For purposes of comparison, the results 
of Edelman and Harsha (ref. 1) are shown and are seen to be relatively 
poor. 

When the complete reaction and eddy breakup models are 
run assuming a confined jet , wherein the axial pressure gradient is 
computed by the code, both tend to lead to too rapid a temperature 
rise (fig. 3); however, in both cases, the predicted maximum tempera- 
ture is only slightly too high. We plan to run this same case 
employing "turbulent" kinetics using the model developed under this 
grant. 

- /oo» 
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Preliminary results for the Northara experiment are in 
figs. 4 and 5. Computer calculations involving this experiment 
will be made in close coordination with HPB personnel. 
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